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PESFACE 

This  paper  presents  the  results  of  a  study  of  imposing  a  three- 
dimensional  perturbation  on  a  "bajotropic  westerly  current.  The 
purpose  of  the  present  investigation  was  to  present  theoretical 
solutions  to  the  following  basic  problems;   (1)  the  determination 
of  the  three  dimensional  struct\ire  of  the  wave^   (2)  the  determination 
of  the  speed  of  propagation  of  the  wave;   (3)  the  value  of  considering 
vertical  motion  in  the  wave  motion. 

Undertaken  as  the  thesis  requirement  for  the  degree  of  I'laster  of 
Science  in  Aerology,  this  peeper  was  prepared  at  the  United  States  Haval 
Postgraduate  School,  Monterey,  California  during  academic  year  1951-1952. 

The  author  is  particularly  indebted  to  Associate  Professor  F,  L» 
Martin  of  the  Department  of  Aerology,  for  his  advice  and  guidance  during 
the  entire  preparation  of  this  study. 
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I.      DTTRODUGTIOKT 

* 

The  theoretical  study  of  the  "behavior  of  large-scale  flow  patterns 
in  the  extratropical  zones  of  the  earth,  associated  with  migratory  long 
waves  in  the  "belt  of  prevailing  v/esterlies,  occupies  an  important  place 
in  meteorological  literature. 

Many  informative  theoretical  and  synoptic  studies  of  the  "behavior 
of  these  large-scale  flow  pa,tterns  have  "been  made  since  V,  Bjerknes  r5  J  » 
in  1916,  advanced  a  theory,  based  on  general  hydrodynamic  considerations, 
that  cyclones  originate  as  dynamically  unsta'ble  wavelike  distur"bances  in 
the  westerly  current.   In  1923 »  H.  H.  Clayton  (^7^   ma,de  the  first 
systematic  use  of  zonal  index  "by  means  of  the  trend  of  mean  pressure 
versus  latitude  to  determine  increasing  and  decrea-sing  zonal  circulation. 
By  1933*  '^^'   Bjerknes  and  collaborators  ^P-TJ  had  computed  the  zonal 
velocity  distri"bution  from  zonal  temperature  and  pressure  distri"bution 
at  the  ground,  showing  the  maximma  horizontal  velocity  normally  occurs 
just  below  the  tropopause  and  concluded  that  there  were  2-3  waves  around 
the  60°  latitude  circle  and  slightly  above  ^  around  the  30*^  circle.   In 
193? »  J.  Bjerknes  r3.28j  offered  a  simple  explanation  for  the  displace- 
ment of  perturbations  superimposed  on  zonal  distribution.  Variations 
in  latitude  and  pressure  perturbations,  with  simultaneous  transport  of 
air,  result  in  changes  of  amplitude  or  longitudinal  displacement  of 
long  waves.   This  result  was  obtained  by  combining  the  equations  of 
motion  with  the  equation  of  continuity.  Rossby  (23j  pointed  out  two 
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"basic  objections,  first,  the  omission  of  acceleration  terms  vjhich 
may  be  of  the  same  order  of  magnitude  as  the  centrifugal  and  de- 
flecting force,  and  second,  the  displacement  of  observed  flow  patterns 
caused  by  isallobaric  systems  resulting  from  divergence  associated  with 
initial  wind  distribution,  Rossby  C27J  extended  the  J,  Bjerknes   paper 
reasoning  tha,t  since  the  depth  of  the  atm.o sphere  was  much  smaller  than 
the  horizontal  scale  of  large  disturbances,  it  was  permissible  to  treat 
the  atmospheric  motion  as  two-dimensional,  non-divergent  flow.  Further, 
he  deduced  that  between  the  longer  westerly  moving  waves,  there  existed 
an  intermediate  wave  length,  such  that  its  perturbation  remained 
stationary.  The  speed  of  propagation  of  moving  waves  is  dependent  on 
the  strength  of  the  westerlies  and  the  wave  length,  and  was  based  on 
the  considera.tion  of  a  change  in  vorticity  in  a  vertical  air  column 
which  has  been  displaced  from  one  latitude  to  another.  He,urwitz  ^16,  if] 
extended  Rossby' s  1939  solution  by  studying  the  effects  of  lateral  limits 
of  the  wave  and  a  variable  Ooriolis  parameter  on  a  curved  earth, 

iJamia.s  and  Glapp  C22I  investigated  Rossby 's  v/ave  speed  formula 
statistically.  Their  verification  was  generally  successful  and  they 
concluded  that  amplitude  and  zonal  velocity  were  both  important  in 
determining  the  stationary  wave  length, 

V.  Bjerknes  had,  by  this  time,  fotmd  the  number  of  waves  around  the 
^5°  latitude  circle  to  be  between  foux  and  six.   In  19^,  J.  Bjerknes 
and  Holmboe  Q^l  derived  a  formula  anala^ous  to  that  of  Rossby  for  the 
more  general  baro tropic  atmosphere  which  agreed  q.ualitatively  with  that 
of  J.  Bjerknes'  of  1937. 
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By  19^7.  Charney  fs]  in  his  'baroclinic  atmosphere  theory,  based 
on  J.  Sjerknes  &jad  Holmboe's  theory  of  wave  motion  in  a  'bajoclinic 
atmosphere,  introduced  a  solution  to  the  barotropic  wave,   ,/ith  this 
theory,  he  obtained  a  solution,  based  on  the  principle  that  a  wave 
will  travel  with  such  a  speed  that  the  pressure  tendencies  8.rise  from 
the  displacement  of  the  pressure  pattern  in  accordance  with  the  field 
of  horizontal  divergence.  Gharney  found  this  physical  explanation 
acceptable,  but  the  dynamics  of  the  wave  could  not  be  solved  by 
analytic  methods  based  on  semi-empirical  considerations  of  gradient 
wind,   Gha.rney  y.9^7j  integrated  the  e«iuations  of  motion,  assuming 
motion  to  be  adiabatic,  and  thus  determined  the  speed  of  propa,gation, 
sta-bility  criteria  and  the  three-dimensional  wave  structure.   In  19^8, 
Slis-ssen  y.2j  generalized  the  Hossby-Holmboe  wave  formula,  but  he 
neglected  vertics,l  acceleration  and  the  Ooriolis  parameter  when  it 
appeared  in  multiplicative  combination  with  the  vertical  velocity. 
In  19^9,  Gha,rney  and  Eliassen  (['7']  extended  numerical  analysis  by  use 
of  the  qxiasi-geostrophic  approxims-tion. 

In  a  field  where  horizontal  divergence  was  so  small  that  the 
absolute  vorticity  of  air  parcels  was  conserved  and  the  scale  of 
motions  was  large  enough  that  the  Ooriolis  parajneter  variation  with 
latitude  was  important  in  determining  wave  dimension,  GressmanL9j 
scrutinized  the  Rossby  wave  velocity  formula  and  found  retrogression 
in  a  series  of  quasi-stationary  long  waves  began  v/hen  the  basic  zonal 
velocity  decreased  in  speed  or  shifted  south.   Such  retrogression  in- 
creased the  long  wave  number  (i.e.,  formed  a  new  cold  trough). 
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Gontinuiiis  his  previous  investigation,  Oharney  ^?,  8j  found  that 
the  mean  flow  with  resx:)ect  to  pressure  could  "be  descrihed  as  two- 
dimensional  non-divergent  in  an  equivalent  harotropic  atmosphere. 
He  used  this  method  as  a  basis  for  partieAly  successful  twenty-four 
hour  forecasts  of  the  500  mh  flow  pattern  "by  numerical  integration 
of  the  non-linear  vorticity  equation.   The  simplified  equations 
derived  expressed  conservation  of  the  vertical  component  of  ahsolute 
potential  vorticity.  The  hydrostatic  equation  together  with  the 
individual  rate  of  change  of  vorticity  evaluated  geostrophically, 
reduced  to  a  single  partial  differential  equation  for  xDressure 
tendency.  The  hasis  of  correspondence  between  the  barotropic  and 
baroclinic  atmosphere  lay  in  the  notion  of  an  equivalent  barotropic 
atmosphere,  in  which  horizontal  motion  at  a  particular  level  approx- 
ima.ted  the  actual  motion. 

It  can  readily  be  seen  that  a  great  amount  of  effort  has  been 
placed  on  the  long  wave  theory  in  an  effort  to  gain  a  forecasting 
tool.  The  use  of  vertical  motions  has  suffered  in  comparison,  but 
primarily  because  of  the  sparsity  of  observations  and  the  necessity 
to  rely  upon  indirect  information.   Investigations  of  the  vertical 
component  of  motion  were  undertaken  as  far  back  as  1911.  when  Y, 
BJerknes  [13}   described  a  method  for  such  computations.   In  1913 » 
3xner  Q.33 suggested  the  high  correlation  between  9  km  pressure  and 
temperatvire  could  be  accounted  for  by  vertical  motions,  Hesselberg, 


W 


in  1915 »  extended  Exner's  argument  "by  showing  that  the  presence  of 
persistent  vertical  motions  of  1-10  centimeters  per  second  would 
account  for  the  oTaserved  magnitude  of  pressure-temperatiire  changes. 
Jeffreys  [^^j  in  his  1922  study  of  vertical  velocities  necessary  to 
maintain  hailstones,  found  vertical  velocities  of  any  magnitude  at 
all  occurred  only  in  eddies  and  could  he  ignored  in  treating  raeaji 
motion.   In  1927,  Haurwitz  ^13]  showed  that  in  an  atmosphere  of 
constant  lapse  rate,  that  the  maximum  vertical  component  of  velocity 
was  in  the  region  of  msociraum  divergence  and  was  of  the  order  of  eight 
centimeters  per  second.  Srunt,  f^J  in  his  Memoirs  0-93^  found  vertical 
velocities  involving  divergence  were  of  the  order  of  a  few  centimeters 
per  second. 

As  observational  information  increased,  quantitative  analysis  of 
divergence  and  vertics-l  currents  became  more  practicable.  Durst  ^Ll^  , 
in  19^0,  found  vertical  velocities  of  a  few  centimeters  per  second; 
Fleagle  \).3l  t    in  19^5*  Graviola  and  Puertes  p.5jl  »  in  19^7.  assumed 
adiabatic  processes  and  computed  the  vertical  component  of  velocity. 
Panofsky  ["237  »  in  19^6,  from  third  degree  polynomials,  found  values 
of  the  same  magnitude  (i.e.,  a  few  centimeters  per  second). 

.iTiile  vertical  motions  of  several  meters  per  second  have  been 
observed,  their  duration  was  short  and  smoothed  out  over  a  twelve- 
hour  period.   Thus,  the  high  valued  measurements  ^oy   theodolites  as 
observed  by  Soring  (iS^  have  not  altered  the  basic  contention  that 
persistent  vertical  motions  axe   of  a  irruch  smaller  order. 
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Charney  L"J»  i^  assigning  magnitudes  to  vertical  motions  in  19^8, 
explained  that  the  failure  of  the  tendency  equs-tion  in  pressure  prog- 
nostication was  due  to  the  discounting  of  vertical  notion,  Charney  t3 J , 
further  determined  that  forecasts  with  the  quasi-geostrophic  assumption 
could  "be  achieved  only  when  the  distrihution  of  vertical  motion  and  sea- 
level  pressure  distribution  v/ere  known.   Thus,  it  is  necessary  to  extend 
numerical  forecasts  to  permit  forecasting  vertical  motion.   In  19^9* 
Bellamy  £23  descrihed  a  method  of  objective  calculation  of  vertical 
velocities.   In  1951*  Phillips  £25]  indicated  that  the  eq.uivalent 
barotropic  atmosphere  was  incapable  of  providing  distribution  of  vertical 
motion  and  sea-level  pressure  adeq.uately  for  forecasting  the  "weather". 

In  view  of  the  interest  in  vertical  velocity  distribution,  it  is  the 
piirpose  of  the  present  investigation  to  present  a  theoretical  solution 
of  a  two-layer  barotropic,  incompressible  atmosphere  with  three-dimensional 
velocity  divergence  zero,  introducing  vertical  motion  &.s   a  perturbation. 
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II.  FOai-IULATIOlI  MB  USSULTS  OF   THE  PHOBLBM 

A  brief  surunary  of  the  principal  subject  matter  of  the  investigation 
unhampered  by  mathematical  detail  is  now  given  to  set  forth  the  procedure 
followed  and  results  obtained. 

Chapter  III  concerns  the  construction  of  a  barotropic  model.  The 
troposphere  is  characterized  by  a  constant  vertiC8,l  Ispse  rate;  the 
stratosphere  is  assumed  to  be  isotherma^l,  with  constant  zonal  wind 
throughout  the  atmosphere.  The  equations  of  motion  are  given  in  ordinary 
cartesian  coordinates  for  an  incompressible  atmosphere.  Horizontal 
velocity  divergence  is  assumed  to  be  compensated  by  vertical  motion. 
Bounde.ry  conditions  are  formulated.   It  is  shown  that  the  model  is 
consistent  with  the  boundary  conditions  and  the  steady  sta-te  notions 
prescribed. 

The  ls.st  section  of  Chapter  III,  treating  the  actual  motion  as  a 
small  perturbation  superimposed  on  an  undisturbed  constant  zonal  current, 
is  now  discussed.   The  velocity  components,  as  well  as  the  pressure 
component,  are  expressed  as  simple  harmonic  perturbations  of  infinite 
lateral  extent  tra.veling  in  a  west-east  direction  at  constant  velocity. 
These  perturbation  components  are  functions  of  height  as  well  as  of 
horizontal  coordinates  and  time.  The  vertical  velocity,  involving  a 
general  expression  of  the  z-coordinate,  is  used  as  a  basis  for  the 
determination  of  all  other  velocity  components.   To  simplify  the 
mathematical  difficulties,  the  meridional  perturbation  amplitude  is 
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made  proportional  to  that  of  the  zonal  perturlDation,   The  zonal 
amplitude,  in  turn,  by  the  prerequisite  of  zero  of  velocity  divergence, 
is  simply  related  to  the  vertical  motion  amplitude  function,  H(z).   It 
is  shown  that  pressure  is  also  a  function  of  the  vertical  motion  ampli- 
tude function. 

Since  the  vertical  amplitude  is  in  terms  of  H(z),  the  amplitudes  of 
the  other  components  are  determinable  from  H(z).  The  determination  of 

the  proportionality  factor  "between  the  meridional  and  latitudinal  per- 

/ 

turbation  velocities  assumes  an  important  role  in  the  solution  of  the 
problem.  Upon  its  determination  rests  the  dependency  of  the  wave 
velocity  and  the  maximam  allowable  meridional  pert\irbation  velocity 
for  a  particular  majcimum  zonal  perturbation  velocity.  Ey  analogy  to 
Rossby's  wave  velocity,  a  restriction  that  the  proportionality  factor 
be  less  than  unity  is  assumed.   It  is  shown  thst  the  proportionality 
factor  is  dependent  ui^on  the  di scrimina-nt  of  an  auxiliary  q_ua-dratic 
equation.  Following  the  boxindary  conditions  this  discriminant  must 
have  a  value  slightly  in  excess  of  unity.   The  discriminant  value 
must  be  determined  from  an  inequs-lity.   The  choice  of  the  value  of 
this  inequality  plays  a  dominant  role  in  any  further  solution  of 
this  problem. 

The  problem  then  becomes  one  of  determining  the,  speed  of  propa- 
gation of  the  wave  and  the  character  of  the  vertical  perturbation 
which  introduces  the  amplitude  as  a  function  of  height.   The  deter- 
mination of  these  are  dependent  upon  wave  length  and  the  parameters 
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characterizing  the  mean  state  of  the  atmosphere,  namely,  the  height 
of  the  tropopause,  the  vertical  velocity  at  the  tropopause,  the 
zonal  wind,  the  vertical  lapse  rate,  the  mean  latitude  and  surface 
density,  pressure,  and  temperature. 

In  Chapter  IV,  the  form  of  the  perturtation  is  set  forth.   The 
pertui'hp.tion  functions  are  solved  in  the  stratosphere,  and  the  tropo- 
si)here  ("by  approximation).   The  proportionality  constant  "between 
meridional  and  zonal  perturbation  velocities  is  determined.   This 
constant  is  of  great  importance  in  the  determination  of  the  wave 
speed  and  stationary  wave  length,  iimplitude  evaluation  is  determined. 
The  complete  solution  of  the  approximate  case  is  shown  diagramatically. 

In  Chapter  Y,  the  tropo spheric  "o^rturhation  functions  are  solved 
v/ithout  simplifying  the  coefficients  of  the  differential  equation 
satisfied  "by  the  Samplitude  function,  H(z). 

Due  to  harotropy,  no  phase  change  is  introduced  into  the  structure 
of  the  v/ave.  Vertical  cells  do  appear  in  the  x-^erturhation  velocity 
structure  moving  with  the  hasic  v;ave  velocity.   The  perturbation 
fields  of  velocity  diminish  v/ith  increasing  height  and  eventually 
approach  zero.  The  wave  in  the  meridional  velocity  field  lags  90° 
behind  the  wave  in  the  pressure  field  and  the  v/ave  in  the  density  field 
is  in  phase  with  the  pressure  wave.  The  zonal  perturbation  is  in  phase 
with  the  pressure  v/ave,  whereas  the  vertical  velocity  perturbation  leads 
the  pressure  wave  by  90°.  These  phase  relationships  are  shown  in 
Figure  3. 
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III.      TE^  I-LaTHEMATICAL  FORl^ULATIOiy 

1,  The  Atmospheric  Model. 

As  an  approximation  to  the  atmosphere  in  middle  latitudes,  a 
model  of  the  undistiirbed  state  is  characterized  as  follov/s;   (a)  the 
undisturbed  v/ind  is  zonal;  the  speed  is  constant  and  independent  of 
the  coordinates  throughout  the  troposphere  and  stratosphere;  (h)  the 
lapse  rate  of  temperature  is  constant  in  the  troposphere  and  zero  in 
the  stratosphere;  (c)  hoth  troposphere  and  stratosphere  are  harotropic, 
incompressihle  and  frictionless  layers. 

2,  Fundamental  Equations, 

It  is  sssuned  for  purposes  of  mathematical  simplicity  that  the 
curvature  of  the  earth  ca.n  be  neglected.  This  $issumption  is  ve,lid  when 
length  of  v/ave  is  small  comp8.red  with  circumference  of  the  zonal  circle 
along  which  the  wave  moves.  This  planar  motion  is  expressed  in  a 
rectangular  system  of  coordinates  x,  y  and  z  with  x  increasing  to  the 
east,  y  northv/ard,  and  z  vertically  upwa.rd.   Then  u,  v,  w  denote 
corresponding  total  velocity  components  (consisting  of  \indisturbed 
velocity  and  perturbation  velocity).  Further,  denote  tota.1  pressure 
by  p,  angular  velocity  of  the  earth  by-TL,  acceleration  of  gravity 
(assumed  constant)  as  g,  geographical  latitude  by  ^   ,  the  z  and  y 
components  of  Ooriolis  acceleration,  f  andy^  respectively.  The 
density  >0  ,  having  been  assumed  incompressible,  is  constajit  and 


(10) 


"'■         equal   to   the  urxperturlDed  value  y^    .      The  iJulerian  equations  of  motion 
in  ordinary  cartesian  coordinates,    fixed  in  the   earth,   hecome 

dt     '      /o  -b^  J  (2) 

dir  /^   9Z  /  '^     '  (3) 

A  fourth  eq-oation  is  obtained  from  the  assumption  that  the  three- 
dimensional  velocity  divergence  is  zero. 


V3  -V  r  O   . 


(4) 


A  relationship  hetween  density  and  height  in  the  troposphere  (with 
the  assumption  that  g  is  constant  v/ith  hei^^.t)  provides  the  fifth 
equation 


^.^L.--3^^     . 


But  the  density-height  relationship  in  the  stratosphere  is 


(5) 


where  h„   is  height  of  tropopause* 


(6) 


(11) 


^  "u.8v^ 
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3.   The  Boundary  Gonditions. 

The  "boundary  conditions  express  the  following  physical  properties 
of  the  raotion:  (a)  the  perturbation  components  are  continuous  across 
the  tropopause,  that  is, 

^P'    A  u.  =*  An?"  -  A  iD-  =  A  "jo  =  AT  -  O.        (7) 

(h)  the  momentum  vanishes  at  the  limit  of  the  atmosphere. 

^«  Derivation  of  the  Perturbation  3qaa-tion. 

The  fundamental  equations  (1-6)  together  with  the  boundary  con- 
ditions (7),  (8)  impose  the  necessary  restrictions  on  the  theoretical 
model,   ..'ith  the  convention  of  a  cs.pital  letter  for  steady  state  value 
and  a  small  letter  for  the  same  quantity  in  the  perturbed  state,  the 
characterization  of  the  steady  state  is  given  by  the  conditions 

\/=  \A/  =  dU  =  "gP  ^  "BlE  c  O   .  U  =  constant.     (9) 

"3T   9t    9\ 

Equation  (2)  implies  that  the  undisturbed  flow  must  satisfy  the 
condition  of  geostrophic  equilibriiim 

o  =   — L^P  ^  f  U  .  (10) 

^Equation  (10)  implies  tha.tyo  and  P  ore  fiinctions  of  y  and  z  only. 
TSquation  (3)  expresses  the  condition  for  hydrostatic  equilibrium  in 
the  mean  state 

o  ^  -_L  a£   -  Q,  (11) 


(12) 


TJU  i^ 


O     =     "vju  <>.        MMtl   «     "^X' 


0 


with  the  omission  ofjQl),   which  can  he  shown  hoth  empirically  and 
theoreticp.lly  to  he  negligihle  in  comparison  to  the  acceleration  of 
gravity. 

By  mea-ns  of  the  equation  of  state,  and  eq.ue,tions  (4)  and  (5)  with 
prescrihed  values  of  U,  Pq  Tq,  f,  "nd  h^,  the  theoretical  steady  state 
model  is  completely  described.   The  boundary  conditions(7)  and  (8)  are 
also  satisfied  by  our  steady  state  model.  Hence  the  model  prescribed 
is  consistent  with  the  fundsjnental  equations  (1-6)  and  the  boimdajy 
conditions  (7),  (8), 
5.   The  linearized  Squations,  i 

The  motion  to  be  investigated  can  be  treated  as  a  small  perturbation 
with  velocity  components  u,  v,  and  w  superimposed  on  an  undisturbed 
constant  zons.l  current  U,  v/ith  the  assumption  of  infinite  lateral 
extent.  Thus  the  motion  is  described  by 

tt  a  u(x,z,t)+ U,  V  =  v(x,z,t),  W  =  v;(x,z,t)  „  (12) 

Similar  expressions  for  the  density  and  pressure  in  the  disturbed  state 
are 

p  =  p(^^,z,t)+P(y,t.),  Jo=   yOCy.z)  .  (13) 

In  deriving  the  perturbation  equations  the  usual  assumptions  are 
ma.de,   (1)  The  total  disturbed  plus  undisturbed  motion  satisfies  the 
hydrodynamic  equa.tions,  as  well  as  the  undisturbed  motion  alone, 
(2)  The  perturbations  are  so  small  that  terms  of  second  order  in  the 


(13) 


'■♦ 


perturbation  quantities  can  be  neglected  with  respect  to  terms  of 
the  first  order  in  the  perturbation  quantities.   Linearizing  the 
"Bulerian  equations  of  motion  (1-3)  snd  simplifying  by  means  of  the 
steady  state  equations  (9-11) i  the  system  (1-3)  becomes: 


9*       5><. 


(15) 


including  the  velocity  divergence. 


7)uc   H-  JS^iT   =  O  .  (17) 


(1^) 
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IV.  THE  BASOTROPIG  PSRTUSBATIOU 

1,   The  Form  of  the  Perturhation  and  Resulting  Conditions, 

The  perturbation  components,  traveling  in  the  x-  direction  at  a 
constant  velocity,  are  assumed  to  be  simple  harmonic  motions  of  infinite 
lateral  extent.  The  amplitiides  of  the  perturbation  are  included  in  a 
general  expression  of  the  z-coordinate,  as  follows 


u  -   sin  m(x-ct)  r(z);  v  «  cos  m(x-ct)  Q.{z) ; 

(18) 
w  -  -m  cos  m(x-ct)  H(z);  p  =  sin  m(x-ct)  G(z). 


To  simplify  the  mathematical  difficulties,  the  assumption  is 
made  that  F(z)  is  proportional  to  0,(2),  D  being  the  constant  of 
proportionality.  Since  u  and  w  of  equations  (18)  are  related  by  the 
condition  of  zero  velocity  divergence,  it  can  be  shown  that 

F(z)  »  H'(z)  =  (^(z)  D.  (19) 

The  equation  for  the  velocity  of  the  wave  is  derived  by  substituting 
the  perturbation  velocities  (18)  into  the  eqiaation  (15) •   iith  the 
equalities  of  (I9)  known,  the  wave  velocity  equation  is 

c  =  U  -  fD  (20) 

m 

Since  U  and  m  are  real,  c  will  be  either  real  or  complex,  depending 

on  D.  The  determination  of  D  will  be  discussed  fully  in  Chapter  IV-2, 


(15) 


(21) 


(22) 


Substituting  the  perturted  velocity,  pressure  and  mess  fields  into 
the  pertixrbation  equations  (l6)  and  (14)  leads  to 

respectively. 

The  differentiation  of  (22)  with  respect  to  z,  leads  to 

By  differentiation  of  density  from  equation  (5)>  we  have 

Equating  (21)  and  (23)  and  replacing  y^  and  ^y<^_  by  their 
expressions  in  (5)  and  (2^)  lef.ds  to  the  following  ordinary 
differential  equation  for  the  amplitude  function  H(z); 

//- o(z) Ha)  -W^^'C^)  -f^AKJ^^  -  Q^O-'^^^c^)=  9   (25) 

where  IC  is  an  abridged  notation  for 
Using  the  further  abridged  notations, 


i^^^. 


(16) 


eq.ijiation  (25) »    the  tropospheric  differential  eq^uation  becomes 

(!-<=< z)  H^z)-  K  H'c^)  +-C^+  BC(-'<4]  Hc^.>  =o  .  (23) 

2.   PerturlDation  Functions  in  the  Stratosphere, 

In  exactly  the  same  manner  as  described  in  the  preceding  section, 
the  stratospheric  differential  equation  follows  "by   using  (6)  instead 
of  (5)  for  the  density  expression.  Thus,  analogously  to  equation  (2^), 
the  stratospheric  density  derivative  is  ^ 

rr  ^yo/A        ,  (29) 


where  ha^^  RT^^/g  is  the  height  of  a  homogeneous  layer  at  str£,tospheric 
temperature  T-,  .v'e,  therefore,  obtain  the  second  order  linear  differ- 
ential equation  v;ith  constant  coefficients 


"'^^-^^feg^.r^']'^'^-'" 


(30) 


The  auxiliary  eoxiation  of   (30 )    is 


(31) 


whence 


q  -  U-   r  -Va  /  / " ■y^gP nih^    .(? 0^^  .     (32) 

The  determination  of  the  value  of  D  is  dependent  upon  the  value 
assigned  to  the  discriminant  of  the  auxiliary  equation  (32).  By 
analogy  with  the  Rossby  wave  velocity 

(17) 


it  would  follow  thet,  for  Rossby  waves, 

The  only  restriction  placed  on  D  at  this  stage  of  the  investigation  * 
is  that  D^  1,  from  which  it  follows  at  once,  that  the  discriminant  of 
(28)  has  a  value  n -  ,  at  latitude  ^5°»  slightly  in  excess  of  uriity. 


nj^>l. 


Solving  the  discriminant  for  D: 


^^■LCSm/,^^ -Cn^  \)~\-C-^^:)i.^  +K  A\=o,  03) 


let 


(3^0 


Solving  equation  (33)  hy  the  qiiadratic  formula  gives 


In  order  to  expand  the  radical  in  (35)  hy  means  of  the  hinonial 
expansion,  the  condition  that 


P*-' j  <  ^l/(rxX-\)  --   2'n'l/f„t-iy 


(36) 


must  hold. 


*  D  7l  with  a  moderately  long  wave  of  6000  km,  would  give  a  value 
of  U-c  >  100  mps,  so  that  ohviously  D  >  1  does  not  apply. 


(18) 


Sxamining  the   two  possible   cases;      (a)   Y^  <  1,      (Id)   X^     1. 
Case   (a):     Y^<1. 
Squation   (3^)  "beoomes 


so  that 


(2whJ  ^(nl-\)Atj 


17?"  <  \-\-ZYr\\i-^     ,  (36') 


%■ 


Jot   L  =  6000  km,  n-j.  <  1.012.  Then  from  the  inequality  of 

case  (a),  the  expansion  of  the  radical  of  (35)fWith  the  choice  of 

the  negative  si^^n, yields 

Case  (l):  Yj^>l 
Squation  (3^)  becomes 


It  follows  that 


where  n^<  1.000^  for  a  wave  length  of  6OOO  km.   The  value  of  n 


(38) 


determined  hy  case  (h)  is  included  in  that  of  case  (a),  so  that  ya:<.l 
will  he  used  in  the  evaluation  of  D.  Moreover,  all  terms  of  (37)  heyond 
the  second  will  he  considered  negligihle.  Thus  at  latitude  ^5°  where 


(19) 


The  evaluation  of  D,  which  for  a  particular  mean  latitude  and  wave 
numher  remains  constant,  is  dependent  on  the  value  of  n-,  selected 
from  the  inequality  (36').   In  this  investigation,  it  was  assumed 
desirahle  to  have  the  ratio  of  the  ahsolute  msjcimom  velocity  of  the 
horizontal  perturbation  to  that  of  the  vertical  perturhation  as  large 
as  consistently  possihle.  This  ratio,  hereafter  called 

M-\!ijmL.\  -  ^^Ci-w.')  .  (39-) 

can  he  shown  to  "be  a  function  of  n  2,  ^^^  1j»  l^'^e  value  of  D  is  com- 
pared in  Tahle  1,  for  various  vs.lues  of  n  2_  s-nd  wave  number  IT. 


H 

^1. 

m  m 

meters'"-^ 

X  10** 

2m  h2_ 

Yz. 

M 

D 

h 

1.0056 

0.888 

0.0113 

0.107 

0.5 

.376 

5 

1.0070 

1.110 

0.01^1 

0.119 

0.5 

.378 

6 

1.008^+- 

1.332 

0.0170 

0,132 

0.5 

.376 

7 

1.0098 

1.55^ 

0.0198 

0.141 

0.5 

.377 

8 

1.0112 

1.776 

0,0226 

0.151 

0.5 

.377 

9 

1.0125 

1.998 

0.025^ 

0.161 

0.5 

.377 

Constant  parameters:  ^  =  45°»  T  =  218°  A 


Tahle  1, 


(20) 


Having  determined  the  value  of  D  q.ualitatively  in  terms  of  m, 
h.  and  n     ,    the  general   solution  of  the   stratospheric  ecius-tion   (30)   is 


/  2.  J 


with  the  abridged  notation 

By  application  of  the  "boundary  condition  (8),yOar-^o   it  follows 

that  f^Cz) ^  O   i^  3-^<i  only  if  the  coefficient  C  is  made  zero. 

The  particulaj:  solution  satisfying  the  "boundary  conditions  is 


(^0) 


\4(z)  :=Ce^-^'-^-^^    . 


(^1) 


In  the  determination  of  the  pressure  perturbation  in  the  stratosphere, 
the  amplitude  function  Gr(z)  is  re'/»ritten  from  equations  (22)  and  (6),  in 
the  abridged  notation,  as 

From  the  solution  of  (^1),  the  pertur"bation  components  of  the 
stratosphere  are  defined.   In  a"bridged  notation,  they  are 

D 

(21) 


I 


Itx  order  to  solve  (^3)  for  Q-^,    suitable  numerical  values  must  "be 
assigned  to  the  constant  parameters  ^  ,  hj,  1!     ,  ^njax*  '^  ^^°-   ^a.  »  /^  * 
The  following  values  are  selected:  ^  =  45°,  h,p  =  10  km,  T^=  218°A, 
Wmax  =  5.5  cm  sec"^,  K  =  5,  n;^_  =  1.00?,  /^  =  0.001121  gm  cm-3. 
With  assignment  of  these  quantities,  numerical  values  may  he  derived 
for  the  parameters  m,  S^,,  M,  where  m  =  1.11  x  10"*^  meter s~l, 
St  =  .0785  km-^,  M  =  -|-.  aquation  (^3)  yields  Gi  =  -5  x  10^  m*"  sec"^, 
3«   An  Approximate  Solution  in  the  Troposphere. 

To  ohtain  the  solution  of  the  harotropic  wave  by  approximating 
the  value  of  (1  -  «,z)   in  the  tropospheric  equation,  (28)  let 

he  a  satisfactory  approximation.   Substituting  this  approximation 
into  equation  (25) i  the  equation  becomes  a  second  order  differential 
equation  with  constant  coefficients,  namely, 

i;  ^"C2^\     -  H  H'Cz)      -h    (A-i-B^X    ^cz)=o,         (45) 

— T—  I  O      / 

r  =Ci-~e)  g//^,      o<  e<  \  .  ,  (^6) 

Then,  if  Ij     ~  ]^1Z  / n    ,  equation  (4-5)  can  then  be  shown  to  be 


^3  K^h^Q:)'--}}        D"-/ 


(22) 


The  general  solution  of  this  equation  is  dependent  upon  the  discriminant 
obtained  from  the  auxiliary  eq.uation  of  equation  (^7).  By  letting  the 
discriminant  equal  >?-  ,  where 

^  ~fTCD'--0  e^CP'-)) 

and  introducing  the  abridged  notation    5g  =  6/2.  ri^  } 

with  D  known  from  the  str8,to spheric  solution  h,  is  determined. 
Then  the  general  solution  in  the  approximate  case  is 

Waj  --  e^^^C^^e"'^''  +  Cje-^jS-^-]  ,        (49) 

4,   Application  of  the  Tropopause  Boundary  Condition. 

The  method  of  evaluating  the  coefficients  Cj^  and  C-  is  the  same 
for  "both  the  approximate  and  actual  case.  The  approximate  case  is  now 
considered. 

With  the  boundary  condition  of  a  discontinuity  of  the  first  order 
at  the  tropopause  the  perturbation  components  must  be  continuous  across 
the  tropopause, 

a^  r  CL3  ^  cxr^  :.  ^  ,  ^  -  -^,  /f-/|  .  (50) 

Then  from  (^9)  the  first  derivative  with  respect  to  2  is 


H'(.)^-S,^'^^   Cc,0..je^^'^%  C3(,^.;e-3VJ_ 


(51) 


(23) 


so  that 


(52) 


From  the  boundary  condition  (50) »  equate  vertical  velocities  8,nd  also 
horizontal  velocities  at  the  tropopause,  thus 

(53) 

By  solving  these  two  equations  sinultaneously  Oj,  and  C^  are  determined, 
5»   The  Complete  Solution  in  the  Approximate  Case. 

>fith  ^2.   ^^-'i  *^3  determined,  the  perturbations  of  the  troposphere 
are  by  approximation  defined  as 


(2^) 


In  order  to  solve  (5^)»  suita"ble  ntunerical  values  must  be 
assigned  to  the  parameters  of  the  troposphere,  in  addition  to  those 
assigned  in  the  stratosphere  (e.g.,  Section  2  of  this  Chapter),   The 
additional  parameters,  having  selected  values,  are 

T«  =  256«5°A,r*  =  5*18°  km"-*-,   v/ith  assignment  of  these  qu8,ntities  it 
can  be  shown  that  the  constants  C^  and  0^  have  the  values  3»59  and 
-4,98  X  10  ,  respectively.  A  diagram  of  the  perturbation  velocities 
is  presented  in  Pigore  1  for  a  graphical  representation  of  the 
solution  determined  in  this  Cha-pter, 


(25) 


im 


V.  TSS   VARIiiBLE  G03i!TIGIi:^TT  GAS3 

1.   The  Solution  of  the  Pert-urhation, 

With  D  determined  and  constant  throughout  the  atmosphere  for 
any  particular  mean  latitude,  the  second  order  differential  equation  (28) 
is  now  solvahle.  By  suitable  manipulation,  descrihed  in  the  Appendix, 
this  equation  (28)  is  a  case  of  the  confluent  hypergeometric  differ- 
ential equation 


{55) 


Equation  {55)    is  satisfied  "by  the  functions 

The  confluent  hypergeometric  equation  then  has  the  general 
solution 


m   =  C:,f,  f  C^%     .  ^^7) 


(27) 


(56) 


To  test  the  convergence  of  f^  -fC^    b  w}   »  compare  to  the 
series  Q^    ,  where  X--*^^  \  K^l  *      It  can  'be   shown  that,  termwise, 
tf^  <  ^^       and  orders  higher  than  the  second  ca-n  be  dropped  with 
accuracy  to  the  first  five  significant  figures.  Then 

where    a/b  =  [/ O^ -\    • 

Testing  the  convergence  of 

after  the  first  12  terms,  the  series  is  more  rapidly  convergent  than 
where  Xs  -  <V^ /^/.   In  the  range  of  values  used  in  this  study,  however, 
there  is  negligihle  error,  if  series  is  shortened  to 

L-y)'^'"'^   ^a  ^  /  +  g-  ^^\      +.  J.  Ca-bfOfa-bf^)/^^      (59) 


2-b     2 


Cp.-^)(3-b) 


The  general  solution  in  the  troposphere  from  equation  {"^1)   is 

Q<-^)'"^  V\^Q^d^±L^   ^j.^^^^/)^jH-2)  J   (60) 


(28) 


2.  Application  of  the  Tropopause  Bo-andary  Conditions. 

After  the  manner  of  evalue.tin^-  coefficients  in  Chapter  IV, 
eq.\aating  the  vertical  and  horizontal  velocity  across  the  tropopause, 
two  resultant  equations  can  he  solved  simultaneously  for  C  and  0-  . 

3.  The  Complete  Solution, 

•iith  Cj^and  C,  determined,   the  perturhations  in  the  troposphere 
are  defined  as 

oT^^    -m  (CiVI    i-  Cj  V^)  Cos  mCK-ct')^ 


I       V>vy\      \  ^ 


(29) 


In  (61)  the  primes  denote  differentiation  with  respect  to  z.  The 
values  of  parameters  to  "be  inserted  in  (6I)  are 


D  =  .377,  m  =  1.1  X  10-8  Q^-l^   oC=  1.9  X  10-"^ 


cm 


-1 


K  =  1.06^  X  10"  cm--'-,  so  that  )p   can  be  evaluated  for  any  value  of  z. 

In  particular,  at  the  tropopause  w  has  the  value  j?  =  -.0176,  Like- 

wise  w   =  .^15  X  10  ,  and  is  constant  throughout  the  troposphere. 

The  complete  solution  includes  the  solution  ahove  the  tropopause 
which  remains  precisely  as  developed  in  Chapter  IV-2. 


(30) 


?I.  Discussion  05"  HSSULTS 

It  would  appear  that  the  variations  of  the  value  of  n2  determined 
"by  the  discriminant  of  the  stratospheric  auxiliary  equation  (32)  plays 
the  most  important  role  in  this  entire  investigation.  By  variation  of 
this  parameter  alone,  the  wave  velocity,  the  stationary  wave  length 
and  the  ratio  of  horizontal  to  the  vertical  maximum  velocity,  can  be 
varied,  ,;ithin  the  allowable  range  of  this  parameter,  wave  velocities 
determined  "by  Rosshy  are  included. 

G-ravity  and  sound  waves  are  also  included  in  these  solutions  for 
appropriate  values  of  n.  The  values  of  n  correspond  to  the  inequality 

ri"^  ^  I   <  Z  tv\  V\     z  ^rrhj,/^        »  with  values  of  the  wave-length,  L, 
suitably  chosen  to  be  representative  values  of  wave-length  for  sound 
and  gravity  waves,  respectively.   It  will  be  remembered  that  Gharney 
^5 J  >  in  dealing  with  the  problem  of  the  unstable  baroclinic  wave, 
first  eliminated  these  short  stable  waves,  which  in  the  present  paper 
is  equivalent  to  restricting  the  condition  to  j  <  1, 

A  further  possible  investigation  is  suggested  by  this  study.  By 
limiting  the  lateral  extent  of  the  disturbance,  larger  values  of  the 
velocity  components  are  consistent  with  the  small  values  of  vertical 
divergence  used  in  this  study,  the  values  of  2i^    and  ^^^ 
would  then  partially  compensate  one  another,  giving  the  proper  (small) 
magnitude  to  ^-  ^  -  -  2idf  . 

Cr  £ 
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APPSilDIX  I. 

The  transformation  of  eq.\iation  (28)  to  the  form  of  a  confluent 
hypergeoraetric  differential  equation  follows, 
aquation  (28)  is  restated  here. 


0-^z)H(z)''ICHCz.)  +(A  +  60-^z)]  Hr^;-o. 


(28) 


Divide  by  (A+3): 


Divide  ^^(f- 8o<2  ^/v  a.  g'\  :^  O         *^®^  multiply  by  (A-f- B) 


(28.a) 


A+B 


H(z) 


Wfe)  t  6^f0)Wa)  =  O    (28.b) 


Momentarily,  consider  the  first  term  only 


{ 


A±£,  -/8fl£,\A±:a. 


8 


B 


( 


A-^6 


3 


(28.c) 


{-^5) 


Expanding, 


I  + 


6o<z 
A +8 


.1' 


(28.d) 


if  it  is  assumed  that  higher  order  terms  negligible. 
Equation  (28, h)  together  with  (28.c,d. )  yield 


(28.e) 


Collecting  terms,  (28. e)  is 


[l-^]H?.)-t^|i+^H^, -^[/'^8]h 


6)  -O. 


(28.f) 


Let 


then 


Afd 

and  equation  (26, f)  is 


A    ^^   I 


5Hc2)  ^  K0^4-!)  Mc^)  4(A4B)H0t)-o,    (28.g) 


(36) 


where 


Squation  (28. g),  in  terms  of  §  ,  becomes 


=  O 


(28.h) 


Dividing  (28. h)  by  [^'V^te]  yields 


^^^)   f 


c.%)(^)^'^)^<^nt^^^^t'-"' 


By  lettlne   k-  =  KBCA+BJ^-^ 


eqtu?.tion  (28,1)  becomes 


y 


B 


(37) 


OThen  if      rif  -  ^, 
equation  (28, j)    is 


1  r*L>)  ^^YA^  -  5)  L'(^)  ^     iitB)!  L^h^)  :.o.       (28.k) 

Divide  equation  (23. k)  by  I'*  and  substitute  the  equality  )fi  s  f  U       ,   and 
the  equation  "becomes 


eq-uation  (23.1)  is  identically  equation  {55)*   which  is  the  form  of 
the  confluent  hypergeometric  differential  equation. 
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